Prepared for submission to JCAP 



Neutron star solutions in perturbative 
quadratic gravity 



Cemsinan Deliduman, afe K. Y. Ek§i c and Vildan Kele$ c 

"Department of Physics, Mimar Sinan Fine Arts University, Bomonti 34380, Istanbul, Turkey 
^National Astronomical Observatory of Japan, 2-21-1 Mitaka, Tokyo 181-8588, Japan 
c Istanbul Technical University, Faculty of Science and Letters, 
Physics Engineering Department, Maslak 34469, Istanbul, Turkey 

E-mail: cemsinan@msgsu.edu.tr, eksi@itu.edu.tr, kelesvi@itu.edu.tr 

Abstract. We study the structure of neutron stars in R + f5W" h ' R^ u gravity model with 
perturbative method. We obtain mass-radius relations for six representative equations of 
state (EoSs). We find that, for |/3| ~ 10 11 cm 2 , the results differ substantially from the 
results of general relativity. Some of the soft EoSs that are excluded within the framework of 
general relativity can be reconciled for certain values of /3 of this order with the 2 solar mass 
neutron star recently observed. For values of /3 greater than a few 10 11 cm 2 we find a new 
solution branch allowing highly massive neutron stars. By referring some recent observational 
constraints on the mass-radius relation we try to constrain the value of f3 for each EoS. The 
associated length scale y/P ~ 10 6 cm is of the order of the the typical radius of neutron stars 
implying that this is the smallest value we could find by using neutron stars as a probe. We 
thus conclude that the true value of /3 is most likely much smaller than 10 11 cm 2 . 
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1 Introduction 

Einstein's theory of gravity, general relativity, enjoyed impressive observational and exper- 
imental support in the past century. Main successes of this theory were explanations it 
brought to phenomena such as precession of the perihelion of Mercury, bending of light and 
the gravitational redshift of light near massive bodies. Although all of these impressive tests 
were done solely in the solar system, the standard model of cosmology assumes the validity of 
general relativity in all scales up to the very large scale structure of the universe. However, as 
the 20th century was ending, data from distant supernovae la [1-3] were interpreted as evi- 
dence of late time acceleration in the expansion rate of the universe. If we continue to assume 
the validity of general relativity in all scales, then the best fit to observational data requires 
us to introduce a non-vanishing positive cosmological constant into the theory. This is the 
simplest way to follow without changing the basic paradigm. Yet there are several theoretical 
problems related with the existence of cosmological constant (see for example [4-7]). One 
of the most important among these problems is the lack of a quantum theoretical method 
to calculate its inferred value from cosmological data [8]. There are several proposals in 
order to avoid the problems of cosmological constant with alternative routes of explanations 
[9, 10]. These proposals can be collected into a few groups: to explain late time accelerated 
expansion one can either modify general relativity by modifying the Einstein-Hilbert action, 
or add new gravitational degrees of freedom other than the metric to the theory of gravity, 
or change how the matter fields and perhaps the cosmological constant gravitates [11, 12]. 

A modified gravity theory, which is proposed to solve late time cosmic acceleration 
problem, should also be able to pass several tests before it can be considered a viable theory 
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of gravity. First of all, in the weak gravity regime, such a theory should be compatible 
with the solar system tests and table-top experiments. In cosmological scales, other than 
producing the late time accelerated expansion, it should be free of gravitational instabilities, 
and obey constraints of the standard model of cosmology. Such a theory is also expected 
to do well in strong gravity regime: for example it should have solutions for neutron stars 
with mass-radius relation inside the current observational bounds. There are alternatives to 
and generalizations of general relativity theory which have the same predictions in the weak- 
field regime as the general relativity, and also provide an explanation of the evolution of the 
universe in the large scale. Thus, the difference between general relativity and alternatives 
might become prominent in the strong-gravity regime [13]. 

One important family of modifications of Einstein-Hilbert action is the f(R) theories of 
gravity (see reviews [14-17] and references therein). In such theories one uses a function of 
curvature scalar as the Lagrangian density. The f(R) term must have a lower order expansion 
in Ricci scalar in order to include general relativity as, perhaps, a weak-field limit of it. Such 
models of gravity can be made to pass Solar System tests, explain the late-time accelerated 
expansion of the universe and also work well in the strong-field regime. In a previous work 
[18] a simple version of f(R) gravity theory, in which f(R) = R + aR 2 , is studied by two 
of the present authors. It is curious that predictions of such theories can be shown to be 
equivalent to scalar-tensor theories of gravity [19] , which have been analyzed throughly since 
the seminal paper of Brans and Dicke [20]. However, if we modify the Einstein-Hilbert term 
by contractions of Ricci and Riemann tensors, R^ v R^ iP and R^ U paR^ upa , then we have an 
alternative gravity theory with independent predictions. 

Inspiration for such an alternative gravity theory may come from string theory. Absence 
of ghosts in low energy string theory in flat backgrounds requires the quadratic corrections 
to Einstein's gravity to be of the Gauss-Bonnet (GB) form [21]: 



However, the Gauss-Bonnet term will not contribute to the equations of motion, because it is 
equivalent to a total derivative in four dimensions. This means that variations of R iW p (J R} LVpa 
can be expressed in terms of variations of R? and Rp U R pv ' . Contraction of Riemann tensors, 
could have relevance only in the cases related to quantum gravity [22, 23]. Therefore for the 
classical physics applications one can take the action of our string-inspired gravity as 



where a and (3 are free-parameters, taken independent of each other. However, special 
relations between a and f3 exists and such cases correspond to some special theories: for 
example, f3 = —3a correspond to the Weyl tensor squared modification of general relativity 
[24, 25], and (3 = —4a has unique energy properties [26]. 

As mentioned above, the absence of ghosts in the low energy string theory in flat back- 
grounds requires the quadratic corrections to Einstein's gravity to be of the form given in 
equation (1.1). This means that the alternative gravity theory defined by (1.2) cannot be free 
of ghosts and this brings up the question of stability of the neutron star solutions discussed 
in this paper. However, all such modified gravity theories, bar the theory defined with (1.1), 
suffer from the same problem, the solution of which is beyond the scope of this paper. Our 
aim in this paper is to see the possibility of neutron star solutions in a specific modified 
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gravity theory defined by (1.2) with a = 0, to analyze solutions with mass-radius relations 
obeying the observational constraints, and to see if there are problems in this gravity model 
other than the well known stability issue. When combined with the results of [18] this analy- 
sis will allow for a comparison of the effects of R 2 and R^R^ terms in (1.2) on the structure 
of neutron stars. 

We further note that the mass dimensions of quadratic corrections to Einstein-Hilbert 
action in (1.2) is order [L] -4 . It is possible to add terms with the same mass dimensions to 
the above action. In fact referring to the quantum gravity arguments in [23], we could also 
add the following dimension [L]~ 4 terms to the action: 

OR and VV^ . (1.3) 

However, as it is pointed out in [27] these terms do not contribute to the field equations and 
therefore in the context of the present paper they are also redundant. 

The effect of the value of a, while (3 = 0, on the mass-radius (M-R) relation of neutron 
stars has been studied, for a polytropic EoS in [28], and for realistic EoS in [18]. The latter 
work constrained the value of a to be |a| < 10 10 cm 2 . In an other study, Santos analyzes 
neutron stars in this model of gravity for a single EoS of ideal neutron gas [29] and for a 
specific choice of ft = —2a, y/a = 0.96 km. In that work, the author shows, for this restricted 
choice of EoS and parameters, that stable configurations of neutron stars are possible even 
for arbitrarily large baryon numbers of the neutron star. 

Our approach in this work is different than [29] in two ways: (i) Rather than choosing 
a fixed specific value for /?, while a = 0, we study its effect, as a free parameter, on the 
M-R relation and constrain its value by referring to observations, (ii) As the interaction 
between nucleons can not be neglected, the EoS of ideal neutron gas can not provide realistic 
M-R relations. We thus employ six different realistic EoSs corresponding to a variety of 
assumptions for the interaction between nucleons rather than the very restrictive EoS of 
ideal neutron gas. 

In the present work, we adopt an alternative theory of gravity, in which the Einstein- 
Hilbert action is modified with the term R^R^ only. In §2, we assume a perturbative form 
of our alternative gravity model and obtain the field equations. The reason of perturbative 
approach is that the equations of motion derived from this alternative gravity theory are 
fourth order differential equations, and their treatment in four dimensions is problematic. To 
obtain the modified Tolman-Oppenheimer-Volkoff (TOV) equations from the field equations 
we also assumed perturbative forms of metric and hydrodynamical functions. In §3, we solve 
the structure of neutron stars in this gravity model for six representative equations of state 
describing the dense matter of neutron stars. We plot the mass-radius relations of neutron 
stars for /3 changing in the range ~ ±10 n cm 2 . This way the value of the perturbative 
parameter f3 is constrained by the recent measurements of the mass-radius relation [30] and 
the observed 2 solar mass neutron star [31]. We identify that (3 ~ ±10 n cm 2 produces results 
that can have observational consequences. Lastly, the results of the numerical study and the 
significance of the scale of the perturbation parameter is discussed in the conclusions. 

2 Modified TOV equations 

The approach that we are going to take in this paper is in the spirit of [18], although not 
in details. The alternative theory that we are going to analyze is not the full quadratic 
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gravity, but part of it whose analysis would complement the analysis reported in [18]. In the 
appropriate units, the defining action of our alternative theory is given as 



S = J d^x^-g (R + PR^BT) + S matter , (2.1) 



where we already set G = 1 and c = 1. This will also be the case in the rest of this section. 
We note that, in this study we are using the metric formalism of gravity in which matter only 
couples to the metric, and the Levi-Civita connection is a function of the metric. Variation 
of the action (2.1) with respect to the metric results in the field equations, 

SirT^ = + p (-^g^RabR ab + V p V p R^ 

+(3 (-VvV^R - 2R^ va R a ° + ^DRg^ . (2.2) 

where G^ v = R^ v — \Rg ilv is the Einstein tensor. The energy-momentum tensor, T^, is the 
energy-momentum tensor of the perfect fluid. 

In the second step we derive the hydrostatic equilibrium equations within the framework 
of the gravity model considered. The hydrostatic equilibrium equations, obtained and solved 
by Tolman-Oppenheimer and Volkoff [32, 33] within the framework of general relativity, are 
commonly called TOV equations. We use the same nomenclature in this paper, although the 
hydrostatic equilibrium equations in this gravity model will turn out to be quite different. 

As in the case of general relativity we choose to work with a spherically symmetric 
diagonal form of the metric and metric functions depend only on the radial coordinate r: 

ds 2 = -e^dt 2 + e 2X ?dr 2 + r 2 {d9 2 + sin 2 Odcf 2 ) (2.3) 

When this metric ansatz and the energy-momentum tensor for perfect fluid substituted 
into the field equations (2.2), the resulting equations will contain metric functions cp(r) and 
A(r), hydrodynamic quantities P{r) and p(r), as well as their derivatives with respect to r. 
The presence of these higher order derivatives precludes expressing the equation in terms of 
hydrodynamic quantities only and consequently deriving the modified TOV equations. In 
order to reach that aim we are going to use the perturbative approach [34, 35] where general 
relativistic solution is taken as the zeroth order solution of the field equations (2.2). As 
mentioned, this method had already been applied to f(R) models of gravity via perturbative 
constraints at cosmological scales [36, 37] and neutron stars with f(R) = R + aR n+1 [18, 28]. 
In the perturbative approach, g^ v can be expanded perturbatively in terms of /3: 

9,w = 9$+P9$+0(P 2 ) (2.4) 
Accordingly, the metric functions must also be expanded in terms of f3 as 

0/3 = + H and A/j = A + H , (2.5) 

and the hydrodynamic quantities of the perfect fluid are expanded perturbatively as: 

pp = p + Ppx + • • • and Pp = P + PP X + • • • . (2.6) 
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In these equations functions without subscript f3 are general relativistic solutions. In the 
present perturbative approach they are "zeroth order" solutions and obtained from Einstein 
equations, 

y± — ±i /\ ji 2 , 

-2 



-8-rrp 
8vrP 



+ e" 



2A (l-2rA>~ 
2A (1 + 2r0>~ 2 , 



(2.7) 
(2.8) 



In order to determine the first and the second TOV equations, we now evaluate tt and 
rr components of field equations (2.2) in terms of metric functions (2.3) and hydrodynamic 
quantities P(r) and p(r). The tt component of the field equations turns out to be 
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and the rr component of the field equations is 



W (P ~ P) 



(2.9) 
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-2A 
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+ (30" - 30' A' + (0') 2 - 2 A') P 



(2.10) 



Here, prime denotes derivative with respect to r. Note that since terms multiplied with (3 
are already first order in /5, metric and hydrodynamic functions that appear in those terms 
are taken as zeroth order general relativistic functions. In all the other terms we have full 
functions. However, when we integrate TOV equations numerically, we are going to keep 
terms only up to first order in (3. 

We now define a mass parameter nip by the relation e~ 2A/3 = 1 — mp/r so that the 
solution for the metric function A / g(r) would have the same form of the exterior solution. 
Here, again we take mp to be expanded in /3 as mp = m + /3m + • • • where m is the general 
relativistic zeroth order solution. It is given in terms of p(r) as 



m 



8ir / p (r) r 2 dr 



Taking the derivative of mass parameter, m^, with respect to r we obtain 



-2X, 



dr 



{2r\' p -l + e 2X e). 



(2.11) 



(2.12) 



In the perturbative approach all terms which are multiplied with f3 in (2.9, 2.10) can be 
rearranged in terms of zeroth order general relativistic expressions by ignoring higher order 

terms. Using the expression (2.12) for — — — in tt field equation (2.9) we obtain the first 

dr 

modified TOV equation as 



dm p 
dr 



8irp/3r 2 — fir 2 



3(l-f)P"-f(pr + 3^-f) pi + \(Pr + ^)p> 



+p 2 + Pp 
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p2 r 3 + m + 2 p m 



(2.13) 
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To derive the second modified TOV equation we use the continuity equation of energy- 
momentum tensor, V^T^ U = 0, which is equivalent to the hydrostatic equilibrium condition, 



dP t 



P 



dr 



ipp+Pp] 



dr 



(2.14) 



We, therefore, need to read — — from the rr field equation (2.10) and use in the above 

dr 

expression to obtain the second modified TOV equation. After some straightforward algebra 
we find it as 



2(r-m /3 ) 



dr 



mr- 



jir 2 



p 2 + Pp-2^(p + P) 

P+P ( p2 r 3 + wP 



2(r~m) 



+ 2Prnj 



(2.15) 



Note that one gets the original TOV equations in general relativity when one sets /3 to 
zero and transforms m — > 2m. The modified TOV equations, (2.13), (2.14) and (2.15), are 
complicated and therefore similar to the case in general relativity they are solved numerically. 
We explain the numerical analysis and present the results of it in the next section. 

3 Solution of the neutron star structure 

In this section we numerically solve the modified TOV equations, namely eqs. (2.13) and 
(2.14), with realistic equations of state (EoSs) appropriate for neutron stars. 



3.1 The Equations Solved 

With the physical constants plugged and m 
become 



2m transformed, Eqns.(2.13) and (2.14) 
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where 



2Gm\ rP" /l 2c 4 3mc 2 \ Gr 2 Pp' 
~)^~ \ 2 + G^P ~ ~^P J 



rc 2 



lp 2 c A lpc 2 2mc 2 pc 2 2mc 2 \2GrP 2 



2 P 2 2 P r 3 P P r 3 P J mc 6 
f 1 pc 2 7mc 2 4c 4 \ GP'r 2 P 
+ V 1+ 2~F~ r 3 P + Gr 2 PJ mc 6 

/ lr 3p mc 2 \ / 2Gm\~ 1 / P\Pp2mG 2 , x 

+ 7^ + ^ + 1 1 H 1 + ^ -r — r- 3 - 4 

The dimension of the coupling constant (3 is length square. Therefore, perturbation is 
actually over the dimensionless quantity (3 = /3//3o> where /3q is the nominal value appropriate 
for neutron stars. It is possible to define it in terms of fundamental constants as follows. The 
typical mass of degenerate stars, including neutron stars, in terms of fundamental constants 
is 

Try 

M = —EL = 1.8482M© (3.5) 

where mpi is the Planck mass and m n is the nucleon mass. Note that this can also be written 
as Mo — a G 3 ^ 2 ?7in where «g — Gm 2 /hc is the "gravitational coupling constant." The radius 
of the neutron star which is about 10-15 km can be given in terms of this mass as three times 
its Schwarzschild radius. As (3q has the dimension of length square, it will be given as the 
square of this typical radius which is the only length scale in the system. Thus 

p = (^^j = 268.128km 2 ~ 2.7 x 10 12 cm 2 . (3.6) 

Accordingly, we can expect that deviations from general relativity (GR), which come through 
the parameter /3, become significant when /3 is a sizeable fraction of This is indeed what 
we observe through the numerical solution of the equations. 

Another way to estimate the order of magnitude of the coupling constant (3 is to see 
that a characteristic value of the Ricci scalar, and the square root of R^R^, is given by 
R ~ GMjc 2 Rl ~ Gp c /c 2 . As f3 ~ 1/R we infer that f3 ~ c 2 /Gp c . In Section 3.6 
we check the validity of the perturbative approach by plotting the dimensionless coupling 
constant /3/A) for a range of static mass configurations. 

In the above equations, the higher derivatives, like P', p' and P" are calculated by using 
the TOV equations obtained in GR. As such terms come only in the perturbative term, the 
error in employing GR instead of the full theory is of the order of f3 2 . 

3.2 Equation of State 

In order to solve the hydrostatic equilibrium equations (3.1) and (3.3) we must supplement 
them with an equation of state (EoS) defining the microscopic physics of neutron stars (NSs). 
The EoS of NS matter at the inner core where most of the mass resides is not well constrained. 
Different EoSs lead to different mass-radius (M-R) relations. As a result we have to solve 
the NS structure for a number of EoSs in order to show that our basic result, that \(3\ < 10 12 
cm 2 , does not depend on the EoSs. In this work we present results for six representative EoSs 
including typical EoS involving only nucleons as well as EoS involving nucleons and exotic 



-7- 



matter. We do not make any analysis on strange quark stars as these are not gravitationally 
bound objects. We use an analytical representation of log/^logP) for all the EoSs obtained 
by fitting the tabulated EoS data following the method described in [38]. These EoSs are 
FPS [39], AP4 [40], SLy [41], MSI [42], MPA1 [43] and GS1 [44]. The physical assumptions 
of these EoSs are described in [45] . 

3.3 Numerical Method 

We employ a Runge-Kutta scheme with fixed step size of Ar = 0.01 km starting from the 
center of the star for a certain value of central pressure, P c . We identify the surface of the 
star as the point where pressure drops to a very small value (10 dyne/cm 2 ) and record the 
mass contained inside this radius as the mass of the star. 

We change the central pressure P c from 3 x 10 33 dyne cm -2 to 9 x 10 36 dyne cm -2 in 
200 logarithmically equal steps to obtain a sequence of equilibrium configurations. We record 
the mass and radius for each central pressure. This allows us to obtain a mass-radius (M-R) 
relation for a certain EoS. We then repeat this procedure for a range of (3 to see the effect of 
the perturbative term we added to the Lagrangian. 

3.4 Observational Constraints on the Mass-Radius Relation 

In order to constrain the value of f3 we have used the recent measurements of mass and 
radius of neutron stars, EXO 1745-248 [47], 4U 1608-52 [48] and 4U 1820-30 [49] (see [46] 
for a description of the method used). We use the 2a confidence contours, shown in all M-R 
plots as the region bounded by the thin black line, on the M-R relation of neutron stars given 
in [30], which is a union of these three constraints.* 

Apart from the above we also use the mass of PSR J1614-2230 with 1.97 ± 0.04 M Q as 
recently measured by [31], as a constraint. It is shown as the horizontal black line with grey 
error-bar. For a viable combination of (3 and EoS, the maximum mass of the neutron star 
must exceed this measured mass. 

These two constraints we employed exclude many of the possible EoSs within the frame- 
work of GR. The value of (3, the coupling parameter of the gravity model studied in the present 
work, is a new degree of freedom which can alter the M-R relation if it takes values of order 
10 11 cm 2 . By using this freedom one can determine the range of (3 values for each EoS that 
is consistent with the observational constraints. 

3.5 The effect of (3 on the M-R relation 

We have determined the M-R relation for each EoSs for a range of f3 values. Results for each 
six representative EoSs are summarized below. In the following we use f3\\ = /3/10 11 cm 2 
instead of (3 as we are mostly interested with f3 of this order. For all EoS we observe that 
the maximum mass of a NS increases with decreasing value of /3n while the radius becomes 
smaller. 

3.5.1 FPS 

The relation between the central density and the mass of the neutron star, p c — M, and 
the M-R relation are shown in Figure 1. For FPS, the maximum mass within GR, is about 
1.8 M which is less than the 2M of PSR J1614-2230 meaning that FPS is excluded within 
GR ((3 = 0). 

*sco [50] for a critic of these constraints . 
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* 1.5 




p c / ( 10 15 g cm" 3 ) 

(a) 




5 10 

R* (km) 



(b) 



Figure 1. The p c — M (panel a) and M-R (panel b) relation for FPS. The solid lines correspond to 
stable configurations for different values of fi. The dashed lines correspond to the unstable configu- 
rations for which dM/dp c < 0. The red line ((3 = 0) stands for the results in GR. The thin black 
line on the right panel shows the 2a confidence contour of the observational measurements of [30]; 
the measured mass M = 1.97 ± 0.04 M Q of PSR J1614-2230 [31] is shown as the horizontal black 
line with grey error bar. The grey shaded region shows where the radius of the NS would be less 
than the Schwarzschild radius. We observe that M max and i? m in increase for decreasing values of (3. 
The comparison with observational constraints on the Figure (see the discussion in §3.5.1) imply that 
—4 < /3n < —2 in this gravity model if FPS is assumed to be the EoS of NSs. 



From Figure 1(b) we find that, for fin = —2, the maximum mass becomes M max ~ 2M . 
Thus FPS can be reconciled with the maximum mass constraint for fin < —2. 

For values of fin < —4 we see that the M-R relation does not pass through the confidence 
contours of the measured [30] mass and radius. This then implies that fin > —4. Together 
with the previous constraint we conclude that FPS is consistent with observations if —4 < 
fill < — 2 in this gravity model. 

Interestingly, we find that for fin < —2 there exists a new branch of stable solution, in 
the sense that dM/dp c > 0, at highest densities which does not have a counterpart in GR. 

For fin = —3 and lower values we find that, for the highest central densities, dm/dr < 
which indicates an unstable star. As a result such solutions are excluded from the Figures 1(a) 
and 1(b). 

3.5.2 SLy 

For SLY the p c — M and M-R relation are shown in Figure 2. For fin > 2 we see that M max 
is less than the measured mass of PSR J1614-2230 [31] and so is not compatible with its 
existence. For fin < —1, the M-R relation does not pass through the confidence contours of 
[30]. These two constraints then imply that 2 > fin > — 1 for the gravity model employed 
here so that SLy is consistent with the observations. 

Again, we see that for fin = —2 (a value which we have already excluded by confronting 
with observations) a new stable solution branch exists for highest densities. For even lower 
values of fi we obtain unstable (dm/dr < 0) solutions for the highest densities. This is why 
we can not extend the solution to fin = —3 for densities greater than 5 x 10 15 g cm -3 . 
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Figure 3. The p c — M (panel a) and M-R relation (panel b) for AP4 EoS. See Figure 1 for the 
notation in the figure. The results are discussed in §3.5.3. 



3.5.3 AP4 

The results for AP4 are shown in Figure 3. It is seen from the M-R relation that for (3n > —2, 
AP4 is consistent with the confidence contours of [30] while the maximum mass becomes lower 
than ~ 2Mq for f3u > 4 and is not compatible with the existence of PSR J1614-2230 [31]. 
Thus we conclude that -2 < /3n < 4 for AP4. 

3.5.4 GS1 

The mass versus the central density, p c — M, and M-R relation for GS1 are shown in Figure 4. 
GS1 is excluded within GR as its maximum mass is well below two solar masses. The 
maximum mass reaches ~ 2Mq for /3n < —3. For /3n < —5, however, the M-R relation does 
not pass through the confidence contours of [30] . The observations than imply the constraint 
-5 < p u < -3 for GS1. 
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For negative values of fin we find that the stability condition, dM/dp c > 0, is satisfied 
for the whole range of central densities considered. As no maximum mass is reached, this 
situation that does not have a counterpart in general relativity. 

3.5.5 MPA1 

The relation between the central density and the mass of the NS, p c — M, and M-R relation 
are shown in Figure 5. The maximum mass of MPAl is above the observed mass of PSR 
J1614-2230 for fin < 12. For fin < 4 the M-R relation does not pass through the confidence 
interval of [30]. This then implies 4 < fin < 12 for MPAl. For larger values of fi the 
maximum mass obtained decreases as is the case for all EoSs. For fin = 8 we see that, apart 
from the usual stable (in the sense that dM/dp c > 0) and unstable branches, there emerges 
a stable branch for the highest central pressures which does not have a counterpart in GR. 
The unstable branch between the two stable branches disappears for even larger values of fi. 
We thus observe that for fin = 12 we can obtain stable configurations with mass exceeding 
2M Q . 

3.5.6 MSI 

The p c — M and M-R relation are shown in Figure 6. This EoS has a maximum mass less 
than 2Mq in GR and does not pass through the confidence contour presented in [30]. In 
order that the M-R relation has a maximum beyond the 2Mq measured mass one needs 
fin < —2 which would take the M-R relation even further from the confidence contour of 
[30]. In order that the M-R relation marginally pass through the confidence contours of [30] 
one needs fin > 2 for which the maximum mass obtained is even smaller. We conclude from 
this analysis that MSI is excluded as a possible EoS for NS in this gravity model assuming 
there is no systematic error in the M-R measurements of [30]. 

3.6 Validity of the Perturbative Approach 

The validity of the perturbative method is justified if the terms that result from the fiR^R^ 
term in the Lagrangian are smaller than the terms arising from R, the GR term. This requires 
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Figure 6. The p c — AI (panel a) and M-R relation (panel b(panel a)) for MSI EoS. See Figure 1 for 
the notation in the figure. The results are discussed in §3.5.6 



that \/3\p c <C c 2 /G be satisfied for each configuration we consider. In Figure 7 we check this 
for EoS AP4 for a range of (3 values. We choose AP4 because it accommodates a large mass 
that could lead to the breakdown of the perturbative approach easily. We see that for large 
masses the dimensionless parameter that measures the perturbative term, \(3\p c G / 'c 2 , reaches 
0.1 which is critical in deciding whether the perturbative approach is justified. This indicates 
that the constraints we provide for f3 may not be accurate for masses approaching 2M for 
\P\~ 10 11 cm 2 for any EoS. 

We observe from the figures in general that variations in the M-R relation comparable 
to employing different EoSs can be obtained for \/3\ ~ 10 11 cm 2 . Using j3 < 10 11 cm 2 gives 
M-R relations that can not be distinguished from the GR results. The range of /3 that is 
consistent with the observations are shown in Figure 8 summarizing the results above. We 
see that |/3n| > 4 is excluded for all EoSs we considered except for MPAl. 
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Figure 7. The dependence of the dimensionless perturbation parameter \(3\p c G/c 2 on mass. The 
result is obtained for the AP4 EoS. We see that the validity of the perturbative approach is question- 
able for large masses and large values of (3. This leads us to conclude that the constraints we provide 
for each EoS may not be as accurate, since the dimensionless perturbation parameter reaching 0.1 
signals the break-down of the validity of the perturbative approach. 



n i i i i 
I 1 FPS 



\ i i i i i i i r 



GS1 



J I I I L 



-| SLY 



AP4 



MPAlJ- 



J I I I I I I I L 



-6-5-4-3-2-1 1 2 3 4 5 6 7 8 9 10 

Pll 



Figure 8. The range of /3n = /3/(10 11 cm 2 ) consistent with the observations. MSI is excluded from 
the plot, because it can not be reconciled with the confidence contours of [30] for the gravity model 
considered though it predicts NS masses well above the measured mass of PSR J1614-2230 [31] for 
/?n < —2. Note that |/3n| > 5 is excluded for all EoSs we considered except for MPAl. The validity 
domain of MPAl extends to very large positive values of (3. We can not include its whole range as 
the perturbative approach applied in this work would break down for such large values. 
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Figure 9. M max changing with j3 for the AP4 EoS. Note that for j3u < —1.8, in this gravity model, 
there is no maximum mass but only stable configurations. 

3.7 Dependence of Maximum Mass on /3 

For all EoSs we observe that the maximum stable mass of a neutron star, M max , and its 
radius at this mass, -Rmiru increases for decreasing values of /3, for the ranges we consider 
in the figures. There is no change in the behavior of M max and R m i n values as /3 changes 
sign. Thus the structure of neutron stars in GR (/3 = 0) does not constitute an extremal 
configuration in terms of M max and R m m. 

In Figure 9 we show the dependence of these quantities on the value of /3 for the EoS 
AP4. We see that /3 = is not special point in figure. 

4 Conclusions 

In this paper we studied the structure of neutron stars (NSs) in a generalized theory of 
gravity motivated by string theory. In the first section we discussed the motivations for 
modifying gravity. In the second section we derived the field equations of this gravity model 
from its action. In the third section we obtained the hydrostatic equilibrium equations in 
spherical symmetry from the field equations by using a perturbative approach in which general 
relativity (GR) stands for the zeroth order gravity model. In the fourth section we solved 
the hydrostatic equilibrium equations for NSs by using numerical methods. In order to solve 
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the equations we have used realistic equations of state (EoSs) that describe the dense matter 
inside NSs and obtained the mass-radius (M-R) relations depending on /3, the free parameter 
of the generalized gravity model considered. These M-R relations are then compared with 
the recent observational measurements of mass and radius of NSs to constrain the value of 
P. 

We have shown that observationally significant changes on the M-R relation are induced 
for j3 ~ 10 11 cm 2 . An order of magnitude smaller values for j3 gives results that are not 
significantly different from what GR predicts. An order of magnitude greater values, on the 
other hand, leads to results that can not be associated with known properties of NSs. For 
such values, the perturbative approach breaks down too. For this selection of EoSs we see 
that none of them are consistent with the observations for /3 < — 5 x 10 11 cm 2 and /3 > 4 x 10 11 
cm 2 . The only exception is MPAl which accepts large positive values like fin ~ 10. We thus 
conclude from this analysis that |/3| <5x 10 11 cm 2 is an upper limit brought by observations 
of NSs [30, 31]. 

We note that, for the the gravity model R + aR 2 , the constraint on a obtained by [18] 
by using NSs confronting the M-R relation with the same observations, is of the order 10 10 
cm 2 , two orders of magnitude smaller than the constraint obtained on (3 in the present work, 
although both a and /3 have the same dimensions ([£] 2 ). This indicates that the R 2 term in 
eq. (1.2) is more effective on the structure of neutron stars than the R^ w R pu term. 

The last observation is very interesting from a theoretical point of view. As mentioned in 
the introduction, theory described by (1.2) is equivalent to Weyl tensor squared modification 
of general relativity for f3 = —3a [24, 25]. That is, from the point of view of Weyl tensor 
squared modification, a and (3 are in the same order. Our perturbative analysis in this paper 
also in some way forces us to make a similar statement: if f3 is taken in the same order as 
of a, then there would not be the problem of the break-down of the perturbative approach. 
What extra we learn from this analysis is that even though the R 2 and R^R^ terms seem 
to be on equal footing in the expansion of a Weyl tensor square term, 

\c^ pu C^ pa = -\r 2 + R^R^ + topological term, (4.1) 

the contributions of them to neutron star physics are not equal. This is pointed out to us 
by the analysis in this paper and it is an important statement just by itself. It would be 
significant if this observation is confirmed and understood analytically. We plan to do this 
in a future publication. 

We find that, some of the EoSs, which do not give M-R relations consistent with the 
observations within the framework of GR, can be reconciled with these observations via the 
free parameter f3 in the generalized gravity model considered in this work. This then brings 
up the question of degeneracy between the EoSs and the free parameter (3. This degeneracy 
does not effect the constraint |/3| < 5 x 10 11 cm 2 which is a bound for all EoSs we considered 
except for MPAl. 

We finally comment that the constraint we obtained is actually the strongest constraint 
we could obtain by using NSs as the experimental apparatus. Another estimate of the 
nominal value (3q of the previous section is as follows: Typical radius of a NS is R* ~ 10 
km, the only length scale in the system. This corresponds to an estimate of curvature 
R pv R^ u ~ R* 2 ~ 10 -12 cm -2 and so the new perturbative term (3R pv R^ u will become of 
order R* and lead to variations on the structure of neutron stars for f3 ~ R 2 ~ 10 12 cm 2 . 
As we obtain such variations in this limit, the value that we should obtain by using NSs, we 
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infer that the actual value of [3 is likely much smaller than this. As we mentioned before, 
deviations from GR are not significant for NSs for values much less than this value. 
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